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Network coding for distributed quantum 
computation over cluster and butterfly networks 

Seiseki Akibue and Mio Murao 


Abstract —To apply network coding for quantum computation, we 

study the distributed implementation of unitary operations over 
all separated input and output nodes of quantum networks. We 
consider networks where quantum communication between nodes 
is restricted to sending a qubit, but classical communication is 
unrestricted. We analyze which IV-qubit unitary operations are 
implementable over cluster networks by investigating transforma¬ 
tions of a given cluster network into quantum circuits. We show 
that any two-qubit unitary operation is implementable over the 
butterfly network and the grail network, which are fundamental 
primitive networks for classical network coding. We also analyze 
probabilistic implementations of unitary operations over cluster 
networks. 

Index Terms —Quantum computing, Network Coding, Quan¬ 
tum entanglement 

I. Introduction 

D ISTRIBUTED quantum computation is computation over 
a network consisting of spatially separated quantum sys¬ 
tems represented by nodes connected by mediating quantum 
systems represented by edges. A serious problem for any kind 
of distributed computation is the bottleneck problem caused 
by the collision of communication pathways between the 
nodes. The bottleneck problem worsens as the network grows. 
Thus it is important to consider how to optimize transmission 
protocols so that the amount of quantum communications is 
reduced. 

In classical network information theory, network coding, 
which incorporates processing at each node in addition to 
routing, provides efficient transmission protocols that can 
resolve the bottleneck problem (TJ. As an example, consider 
a communication task over the butterfly network and the grail 
network presented in Fig. Q] that aims to transmit single bits 
x and y from i\ to 02 and 12 to <>\ simultaneously via nodes 
rii, 712, 7 i 3 and 714. The directed edges denote transmission 
channels with 1 -bit capacity. One of the channels in each 
network (the channel from m to 712 for the butterfly network 
and either the channel from rii to 712 or the channel from 773 
to 714 for the grail network) exhibits the bottleneck without 
network coding shown in Fig. [T| 

Quantum communication with quantum network coding has 
been studied by analogy to classical network coding 0 , 
0 , 0 , 0 , 0 , 0 . fc-pair quantum communication over a 
network is a unicast communication task to faithfully transmit 
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Fig. 1. Network coding for a classical communication task over i) the butterfly 
network and ii) the grail network. Two bits of information x,y E {0,1} are 
given at the input nodes i\ and 22 , respectively, x © y denotes addition of x 
and y modulus 2. 


a fc-qubit state given at distinct input nodes {*1,72, ••• ,ik} 
to distinct output nodes (oi,02,--- ,Ok } through a given 
network. Two examples of 2 -pair quantum communication 
over a butterfly network and a grail network are shown in 

Fig- 0 

In quantum mechanics, the no-cloning theorem forbids 
the creation of a perfect copy of an unknown state. Thus 
perfect multicast communication of an unknown input state 
is impossible. As copying states is a key element of classical 
network coding, classical network coding cannot be simply 
extended to fc-pair quantum communication over the net¬ 
works. Indeed, in the setting where each edge can be used 
for either 1 -bit classical communication or 1 -qubit quantum 
communication, perfect quantum 2 -pair communication over 
the butterfly network has been shown to be impossible 0 , 
0 - However, it has been shown that if each edge can be used 
for either 2 -bit classical communication or 1 -qubit quantum 
communication, perfect quantum 2 -pair communication over 
the butterfly network is possible, if and only if input nodes 
share two Bell pairs 0 . 

Further, if each edge has 1 -qubit channel capacity and 
classical communication is freely allowed between any nodes, 
it has been shown that there exists a quantum network cod¬ 
ing protocol to achieve the 2 -pair quantum communication 
over the butterfly and grail networks perfectly 0 , 0 , 0 - 
This setting is justified in practical situations, where classical 
communication is much easier to implement than quantum 
communication. Moreover in this setting, the correspondence 
between classical and quantum network coding has been 
discovered. Namely, it has been shown that fc-pair quantum 
communication is possible over a network if the corresponding 
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Fig. 2. i) The butterfly network and ii) the grail network with the input nodes 
( i\ and 12 ), output nodes {o\ and 02 ) and the repeater nodes ( 711 , 712,773 
and 714 ). The directed edges E 2 , • • • , Eg represent quantum or classical 
channels. Quantum channels have 1-qubit capacity. There are several settings 
of classical channel capacity for each edge. Our task is to transmit a given 
two-qubit state |input)j l 5 j 2 from i\ to 02 and from 12 to 01 simultaneously 
by using the channels and local quantum operations at each nodes under the 
setting where classical communication is freely allowed between edges. 

/.-pair classical communication is possible over the network 
using linear classical network coding schemes 0, © or even 
using nonlinear schemes 0 . A connection between linear 
classical network coding schemes and measurement-based 
quantum computation has been investigated in 0. However, 
it has been an open problem to determine the possibility 
of /.'-pair quantum communication over networks where the 
corresponding /c-pair classical communication is impossible. 

In /.--pair quantum communication, the output state 
|output) 0i ... 0fc can be regarded as a state obtained by per¬ 
forming a /e-qubit unitary operation U on the input state 
linput)^...^ 

|output) 0i ... 0fe = U |input) il ... ifc , (1) 

where U is a permutation operation. We do not need to restrict 
the fc-qubit unitary operation U in bq.(|T]i to be a permutation 
operation, it can be a general quantum operation. This leads to 
the idea of network coding for quantum computation , which 
aims to perform a quantum operation on a state given at 
distinct input nodes and to faithfully transmit the resulting state 
to the distinct output nodes efficiently over the network at the 
same time. By computing and communicating simultaneously, 
quantum computation over the network may reduce commu¬ 
nication resources in the distributed quantum computation 
scenario. The study of network coding for computation is still 
in its infancy for classical and quantum cases. Network coding 
for classical computation is considered in Col and network 
coding for quantum computation over the butterfly network is 
considered in 0, both in 2011. 

In this paper, we investigate a cluster network, which is 
a special class of network with k input nodes and k output 
nodes, as a first step to apply network coding for more 
general quantum computation. The cluster network contains 
the grail network as its special case. We focus on the setting 
considered in 0, ©, 0, where classical communication is 
freely allowed between any two nodes. We identify the class 
of unitary operations that can be implemented over cluster 
networks in this setting by investigating transformations of 
cluster networks into quantum circuits implementable by using 


quantum communication resources between nodes specified 
by the cluster network. The transformation method of cluster 
networks is also applicable to the butterfly network. It provides 
constructions of quantum network coding for implementing 
any two-qubit unitary operations over the grail and butterfly 
networks, which are the fundamental primitive networks for 
classical network coding. We also analyze probabilistic im¬ 
plementation of V-qubit unitary operations over the cluster 
network to understand the properties of quantum network 
coding for quantum computation when the requirement of 
deterministic implementations are relaxed but that of exact 
implementations are kept. 

The rest of this paper is organized as follows. In Section 

II, we review necessary background information. In Section 

III, we define the cluster network and the implementability of 
a unitary operation over a quantum network. In Section IV, 
we present a method to convert a given cluster network into 
quantum circuits describing unitary operations that are imple¬ 
mentable over the network. In Section V, we show that any 
two-qubit unitary operations is implementable over butterfly 
and grail networks. In Section VI, we investigate the condition 
for unitary operations to be implementable over a given cluster 
network and show that our conversion method presented in 
Section IV gives all implementable unitary operations over 
the cluster networks with 2 and 3 input nodes. Since the 
condition for unitary operations presented in Section VI is 
not based on classical network coding schemes but based 
on properties of quantum operations, it can be used for an¬ 
alyzing /.'-pair quantum communication over networks where 
corresponding fc-pair classical communication is impossible. 
In Section VII, we investigate probabilistic implementation 
of unitary operations over the cluster network and show that 
it is impossible to achieve 2-pair quantum communication 
even probabilistically over a square shaped four-node quantum 
network, where corresponding 2-pair classical communication 
is impossible. A conclusion is given in Section VIII. 

II. Preliminaries 

A. Notations 

The following notation will be used throughout this paper. 


a 

The complex conjugate of a. 

a T 

The transpose of a. 

a t 

The conjugate transpose of a. 

L («) 

The set of linear operators 
acting on the Hilbert space TL. 

Ia 

The identity operator on 1~La- 

U(«) 

The set of unitary operators. 

u c 

The set of unitary operators 
locally unitarily equivalent to 
a two-qubit controlled unitary operation. 

tr 

The trace of a linear operator. 

L {Ha :Hb) 

The set of linear operators. 

L : I-La —>• 'Hb- 

U (H A :Hb) 

The set of isometry operators. 

Sch#£(|V>» 

The Schmidt rank of E 'H.a 

Op #s(M) 

The operator Schmidt rank of 

KC#(U) 

The Kraus-Cirac number of 
a two qubit unitary operator. 
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B. The Schmidt decomposition and rank 

For any vector | i/j)ab in a Hilbert space TLa ®TLb, there 
exist a set of orthonormal vectors {|*}a £ P-A^i and a set of 
orthonormal vectors {\i)b £ Pb}i such that 

\4>)AB = ^2 ^i\i)A\i) b, (2) 

i 

where {A^}i are non-negative real numbers referred to as 
Schmidt coefficients. The decomposition of a vector \iP)ab 
given in the form of Eq. Q is referred to as a Schmidt de¬ 
composition of | iP)ab- Each Schmidt coefficient is equivalent 
to the square root of an eigenvalue of the reduced density 
operator pA on TLa of | i/j)ab given by 

p A = tr B (\ip}(i>\AB), ( 3 ) 

where tr b(X) = JA (bi\sX\bi)B is a partial trace of I £ 
L(TLa ® TLb ) and {| b/)B £ Pb} is an orthonormal basis. 
The Schmidt decomposition is uniquely determined up to 
arbitrary choices of the orthonormal vectors in the subspaces 
corresponding to degenerate Schmidt coefficients. 

The number of non-zero coefficients |{A* > 0}| is called as 
the Schmidt rank of \iP)ab ■ The Schmidt rank of \iP)ab is 
denoted by Sch#^(| i/j)ab) = |{> 0}| in this paper. For a 
vector \i/i)abc in a multiple Hilbert space TLa ® TLb 0 TLc, 
the Schmidt decomposition of \i/j)abc in terms of a bipartite 
devision between the Hilbert spaces Pa®Pb and TLc can be 
similarly defined by introducing a set of orthonormal vectors in 
each devision and the corresponding Schmidt rank of \i/j)abc 
is denoted by Sch#^ b (| iJj)abc)- 

C. The operator Schmidt decomposition and rank 

The operator Schmidt decomposition can be applied to any 
linear operators acting on a Hilbert space TLa 0 TLb ■ The set 
of linear operators L (TLx) forms a Hilbert space with respect 
to the inner product ( M,N ) = dim ^ A .) tr(M^N). Thus we 
can apply the Schmidt decomposition to operators, such that 
for any linear operator M G L (Pa®Pb), there exists a set of 
orthonormal operators {Pi £ L(Hi)}; and {Qi £ L (TLb)}i 
satisfying 

M = J2 X * p i ® Qi’ ( 4 ) 

i 

where {A^} are non-negative real numbers referred to as 
operator Schmidt coefficients l26l of M. The decomposition of 
a linear operator M given in the form of Eq. (Q is referred to 
as the operator Schmidt decomposition of M. The number of 
non-zero coefficients | {A^ > 0} | is referred to as the operator 
Schmidt rank of M. In this paper, we denote the operator 
Schmidt rank of M by Op #g(M). 

D. The Kraus-Cirac decomposition and rank 

A general two-qubit unitary operation U £ U (TLa 0 TLb) 
where TLa = Tin = C 2 can be decomposed into a canonical 
form called the Kraus-Cirac decomposition introduced in lITTfl . 
Oil, Ell given by 

U = (U 0 u >) e i^X9X+yY9Y+zZ9Z) ^ g ^ ( 5 ) 


where u, v!, w and it/ are single-qubit unitary operations and 
X, Y and Z are the Pauli operators on C 2 and y. z £ R. 
In particular, the two-qubit global unitary part U g i 0 b a i{x, y , z) 
of U is defined by 

U g i 0 bal{x,y,z) := (6 ) 

In Eq. the parameters x, y, z in 0 < x < Tr/2 (or 0 < 
x < 7r/4 if z = 0), 0 < y < min{a;, 7r/2 — x} and 0 < z < y 
cover all two-qubit global unitary operations up to the local 
unitarily equivalence (the Weyl chamber 1 13| |). 

The Kraus-Cirac number of a two-qubit unitary operation 
U is defined as the number of non-zero parameters x,y,z 
in U globalise , y , z) and is denoted by KC#(U) in this paper. 
KC#(C/) characterizes nonlocal properties (globalness) of U 
m. The following is the list of classifications of two-qubit 
unitary operations: 

• U with KC#(U ) = 0 is a product of local unitary 
operations and satisfies O P# f, ( U ) = 1. 

• U with KC#(U) = 1 is locally unitarily equivalent to a 

controlled unitary operation and satisfies = 2. 

• U with KC#(U) = 2 is locally unitarily equivalent to 
a special class of two-qubit unitary operations called a 
matchgate ITBl , fl6l . fTTI and satisfies OP#g(U) = 4 

• The rest of two-qubit unitary operations including a 

SWAP operation have KC#(Z7) = 3 and satisfies 

Op#b(U) = 4. 


III. Cluster networks 

We denote the Hilbert space of a set of qubits specified by 
a set Q by TLq and the Hilbert space corresponding to a qubit 
Qi specified by an index k by 'TLo k = C 2 . In our setting 
where quantum communications are restricted but classical 
communications are unrestricted, quantum communication of 
a qubit state between two nodes is replaced by quantum 
teleportation El between two nodes. Since any direction of 
classical communications is allowed, quantum communication 
of a single qubit state can be achieved by sharing a maximally 
entangled two-qubit state between the nodes and the direction 
of quantum communication is not limited. Thus quantum 
network coding is equivalent to perform local operations (at 
each nodes) and classical communication (LOCC) assisted by 
the resource state that consists of a set of maximally entangled 
two-qubit states (the Bell pairs) |4> + ) = ^(|00) + |11)) shared 
between the nodes connected by edges. 

We investigate which unitary operations are implementable 
by LOCC assisted by the resource state for a given network 
where nodes are represented by a two-dimensional lattice. We 
consider that a node represented by Vij is on the coordinate 
of the two-dimensional lattice (i,j) and edges connect nearest 
neighbor nodes. We call these networks cluster networks. We 
first give a formal definition of a cluster network. 

Definition 1. A network G = {V,S,T 1 0} is a (fc, N)-cluster 
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network if and only if, 

V = {i1 < * < k, 1 < j < N} 

X = {vi t i; 1 < i < k} 

O = {vi,N\ 1 <i<k} 

£ = S U K, (7) 



where 


S = l<i<k-l,l<j< N)}, 

K- = {( v i,ji v i,j+1 ); 1 <i <k, l<j<N- 1)}, (8) 

k > 1 and N > 1. V represents the set of all nodes, X and 
O represent k input nodes and k output nodes, respectively. 
£ represents the set of all edges and S and K, represent the 
sets of vertical and horizontal edges, respectively. 


Next we define the resource state corresponding to the 
( k, ./V)-cluster network. We introduce qubits ,5 ',? ;j at node v t ,j 
and Sf +1 j at node Vi+ij to represent a Bell pair correspond¬ 
ing to an edge (vij, fi+ij) G S. Similarly, we introduce 
qubits specified by Kfj at node v l:1 and Kfj +1 at node 
Vij +1 to represent a Bell pair corresponding to an edge 
K,. We denote the set of all qubits in the 
resource state by 1Z = {Sjj,Sf +1 j |1 < i < k — 1,1 < j < 
N}U {Klj, Kfj +1 \l < i < k, 1 < j < N — 1}. The resource 
state corresponding to a cluster network is defined by 

the following. 


Definition 2. For a given (k, N)-cluster network, the resource 
state | < 1?)7?, G FLn is defined by 


= 


i 




J=1 \® + )sij,S} +1 ' 
=1 


(9) 


For example, the (3,3)-cluster network and the correspond¬ 
ing resource state are shown in Fig. [3] Note that the resource 
state for a cluster network represented by Eq. @ is related 
to but different from the cluster states used in measurement- 
based quantum computation ED. While we can convert the 
resource state for a cluster network into a cluster state by 
applying a projective measurement on all qubits at each 
node, a cluster state cannot be converted to the resource 
state for the corresponding cluster network by LOCC. The 
resource states for a cluster network is also closely related to 
a valence bond solid state |[20l introduced in condensed matter 
physics through the projected entangled pair states (PEPS) EE) 
representation for the valence bond solid states. The resource 
states for cluster networks are equivalent to a special type of 
resource states consisting of Bell pairs used for representing 
2D (square lattice) PEPS. PEPS can be represented as states 
probabilistically obtainable by independently performing a 
linear transformation on each node on a resource state. In 
contrast, conditional operations at each node depending on 
the outcomes of measurements in other nodes are performed 
in our network coding scheme for quantum computation. 

Finally we define the implementability of a unitary oper¬ 
ation over a A:-pair network. In addition to resource qubits 
1Z, we introduce input qubits /, at the input node j G X, 
output qubits Oi at the output node v 1: .m £ O, a set of input 


ii) 


input 

qubit 


output 

qubit 



Fig. 3. i) The (3,3)-cluster network with the input nodes X = 
{1^1,1, ^2,1, ^3,1}, output nodes O = {^1,3, ^2,3, ^3,3} and 3 repeater nodes 
{^1,2, ^2,2, ^3,2}, and ii) the corresponding resource state. Note that the 
resource states of the cluster networks are different from the cluster states 
used in measurement-based quantum computation d 


qubits Xq = {/j|l < i < fc} and a set of output qubits 
Oq = {Oi|l < i < k} for a (k, iV)-cluster network. Note 
that each input and output node has only one input or output 
qubit since we concentrate on the implementability of a unitary 
operation, and the state of output qubits is initially set to be 
in |0) G Uo Q - 

Definition 3. For a (fc, N)-cluster network specified by G = 
{V,£,X,0}, a unitary operation U G U fHx Q '■ T~Lo Q ) is 
deterministically implementable over the network if and only 
if there exists a LOCC map F such that for any pure state 
\4>) 6 FLiq, 

r(|'0X'0l ® |$}<$k) = (io) 

where LOCC map F consists of local operations on each node 
and classical communications and U (fH.x Q '■ XLo Q ) is the set 
of unitary operations from FLx Q to FLoq- 

Note that the main difference between this network com¬ 
putation model for implementing a unitary operation over 
a cluster network and standard measurement-based quantum 
computation is that any operations inside each node are 
allowed including adding arbitrary ancilla states in this model 
whereas only projective measurements on the cluster state in 
each node are allowed in measurement-based quantum com¬ 
putation. Thus the full set of implementable unitary operations 
over a (k, 7V)-cluster network is larger than or equal to a set 
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of operations implementable by measurement-based quantum 
computation using the corresponding cluster states converted 
from the resource state for the (k, iV)-cluster network by 
LOCC. 

IV. Conversion of a cluster network into 

QUANTUM CIRCUITS 

We propose a method to convert a (k, A')-cluster network 
into quantum circuits representing a class of unitary operations 
implementable by LOCC assisted by the resource state corre¬ 
sponding to a given cluster network. By using the converted 
circuit, it is easier to construct a network coding protocol since 
a set of implementable unitary operations are represented by 
a set of parameters of the converted circuit instead of a com¬ 
plicated LOCC protocol. The class of implementable unitary 
operations represented by the converted circuit is a subset of 
that over the cluster network in general since this particular 
conversion method does not guarantee to give all possible 
constructions. However, in some cases, the constructions given 
by the conversion methods cover all possible implementable 
unitary operations as will be shown in Section V. 

We define a set of vertically aligned nodes VJ := {}f =1 
and a set of vertically aligned edges Sj := {{vjj, Ui+i,j)}fcTi 
where 1 < j < N. We also define a set of horizontally aligned 
nodes Vj := {vjj }jL\ and a set of horizontally aligned edges 
Ki := {(v i , j ,Vi !j+1 )}fs 1 1 where 1 < i < k. We consider that 
the Bell pairs given for a set of vertically aligned edges Sj 
are used for implementing global unitary operations between 
nodes whereas each Bell pair given for a set of horizontal 
aligned edges ICi is used for teleporting a qubit state from 
node Vj j to node Vi.j+ 1 . 

We show that three types of unitary operations, a two-qubit 
controlled unitary operation, a three-qubit fully controlled 
unitary operation and a single qubit unitary operation, are 
implementable between the nodes in VJ if only one Bell pair is 
given for each edge and LOCC between the nodes is allowed. 
Details of a LOCC protocol implementing three-qubit fully 
controlled unitary operations are given in the next subsection. 
A LOCC protocol implementing two-qubit controlled unitary 
operations has been proposed by l22l and also obtained by 
simply applying the protocol implementing three-qubit fully 
controlled unitary operations, as the two-qubit control unitary 
operations are special cases of three-qubit fully controlled 
unitary operations. 

• A two-qubit controlled unitary operation: A two-qubit 
controlled unitary operation is defined by 

l 

'■='^2\a)(a\i®ul?\ (11) 

a=0 

where l represents the vertical coordinate of the node vij 
of a control qubit and n represents the vertical coordinate 
of the node v n3 of a target qubit, and uif 1 (a = 0,1) 
are single qubit unitary operations on the target qubit. 
If n ^ / ±1, all Bell pairs represented by edges 
between l and n are consumed to implement the two- 
qubit controlled unitary operation. When we do not 
specify the single qubit unitary operations {u^} on the 


target qubit we denote a two-qubit controlled unitary 
operation simply by C; ;ra . In particular, if Un ) = I „ and 
= X, Ci- t n({u^} a - o,i) is called as a controlled- 
NOT operation. 

• A three-qubit fully controlled unitary operation: In addi¬ 
tion to the two-qubit control unitary operations, we can 
perform three-qubit fully controlled unitary operations 
defined by 

l l 

Ci,mA{Un b) }a,b=o,i) ~ ]T X] \ab) (ab\ lm ® , (12) 

a—0 b —0 

where l and m represent the vertical coordinates of the 
nodes vij and v m j of two control qubits, respectively, 
and n represents the vertical coordinates of the node v n j 
of a target qubit, and ui? b \a,b = 0,1) represents single 
qubit operations on the target qubit. (See the next subsec¬ 
tion for details of the LOCC protocol implementing three- 
qubit fully controlled unitary operations.) Note that the 
indices l, n and m should be taken such that l < n < m 
or m < n < l. Similarly to the case of a two-qubit 
controlled unitary operation, we denote a three-qubit fully 
controlled unitary operation by Ci tTn;n when we do not 
specify the target single qubit operations. On the other 
hand, any four-qubit fully controlled unitary, where three 
of the four qubits are control qubits and the rest of the 
qubit is a target qubit, is not implementable on qubits 
that are all in different nodes of VJ in a (k, A')-cluster 
network, if a single Bell pair is given for each edge in 

Sj- 

• A single qubit unitary operation: Obviously any single 
qubit unitary operations can be implemented on any qubit. 

Note that a general three-qubit fully controlled unitary 
operation is not implementable by using a sequence of two 
two-qubit controlled unitary operations that is implementable 
by using vertically aligned Bell pairs in general. This indicates 
that the use of three-qubit fully controlled unitary operations 
enhances the implementability of converted circuits. A three- 
qubit fully controlled unitary operation plays an essential role 
in our network coding protocol over the butterfly network as 
shown in the next section. 


A. A LOCC protocol for implementing three-qubit fully con¬ 
trolled unitary operations 

We show a construction of a protocol to implement a 
three-qubit fully controlled unitary operation on qubits 

located at a set of vertically aligned nodes VJ over the (k. N)- 
cluster networks, where l and m represent two control qubits 
at nodes vij and v m j respectively, and n represents a target 
qubit at node v n j. We present a LOCC protocol to implement 
Cz,m;n assisted by the resource states consisting of the Bell 
pairs corresponding to the vertical edges Sj of the (k. N')- 
cluster networks. 

We consider to implement Ci }m;n on a state of three qubits 
indexed by Qi, Q m and Q n at node vij, v m j and v n j, 
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respectively, and its explicit form is given by 

1 1 

Cl,mA{^n b) }) ■■= E E \ ab )( ab \lm ® Ui ab) (13) 

a—06—0 

where {|ah)} a ,{>={o,i} is the two-qubit computational basis of 
TLq, <S> T~LQ m of the two controlled qubits and u'n^ acts on 
%Q n of the target qubit. 

To show how our LOCC protocol works, we consider an 
arbitrary state of the control qubits by ^ \ a b\ab)im £ Hq, <g> 
'H.Q m where {A„ ( ,} is a set of arbitrary complex coefficients 
satisfying the normalization condition ]P a b X a h \ 2 = 1 and we 
represent an arbitrary state of the target qubit by 1<j>) £ 'Hq„. 
In the following, we show that our protocol transforms the 
joint state of controlled qubits and a target qubit to 

C/,m;n ^ \o,b)l m | (jt) n ~ \ ^ | rt - 

a,6 a,6 

The protocol for implementing three qubit fully controlled 
unitary operations (see Fig. II is specified as follows: 

1) Ancillary qubits indexed by Qy, Q m > are introduced at 
nodes vij and v m j respectively. Set both of the ancillary 
qubits to be in |0). Each of the two states of control 
qubits Qi and Q rn is transformed to a two-qubit state 
by applying a controlled-NOT operation on the control 
qubit and the ancillary qubit at the same node, namely 
applying controlled-NOT operations on Qi and Qy and 
also Q m and Then the joint state of five qubits Qy 
Qv, Qm, Qm' and Q n is given by 

'^\ab\ab)i m \ab)y m '\4>)n- (14) 

a,6 

2) By consuming the Bell pairs corresponding to the ver¬ 
tical edges Sj between vij and v n j and also between 
v m j and perform quantum teleportation to transmit 
the states of qubits Qy and Q m > from nodes vij and 
v rn j to v n> j. A circuit representation of the protocol 
of quantum teleportation represented by T in Fig. Q] is 
given by Fig. [5] We denote indices of two qubits at node 
v n j representing the teleported states from Qy and Q rn > 
by Qy and Qm", respectively. 

3) At node v n j, perform C^ m - n on 'Hq 1 „ 

Then we obtain the state given by 

y: \ab\ab)lm\ab)v'rn"U ( £ b) \<t>)n ■ (15) 

a,6 

4) At node v n j, we apply the Hadamard operations and 
perform projective measurements in the computational 
basis on both 'Hq 1 „ and Hq „ ■ The measurement 
outcomes of qubits Qy and Q m " are sent to nodes vy :j 
and v m j, respectively, by classical communication. At 
each of nodes vij and v rn .,j, if the measurement outcome 
is 0, do nothing, and if the outcome is 1, perform Z 
for a correction on qubit Qi or Q m ■ By straightforward 
calculation, we obtain the state of three qubits Qy Q m 
and Q n at nodes vij,v m> j and v n j, respectively, given 
by 

'Y^Kb\ab)i m u^ b) \4>) n - (16) 

a,6 


Vi,j 


Vn,j 


Vm,j 



Fig. 4. A quantum circuit representation of the LOCC protocol implementing 
a three-qubit fully unitary operation Ci m-n, where the qubits in the first 
shaded region are at the node vij, those in the second shaded region are at 
the node v n j and those in the third shaded region are at the node Vmj. The 
protocol consists of entangling ancillary qubits Qy and Q m / at the nodes V[ j 
and Vmj, respectively, by performing controlled-NOT operations at vi j and 
v m,j, teleporting ancillary qubit states from the nodes vi j and v m j to the 
node v n j represented by qubits Qyi and Q m // by applying a teleportation 
protocol denoted by 7”, applying C\ m . n on controlled qubits Qy/ and Q m // 
and a target qubit Q n at the node v n j, performing Hadamard operations and 
measurements in the computational basis on Qy/ and Q m // at node v n ,j 
and finally applying conditional Z operations depending on the measurement 
outcome on two control qubits Qy Qm at nodes vi j and v m ,j , respectively. 



Fig. 5. A quantum circuit representation of the quantum teleportation protocol 
T, where represents the measurement projected in the Bell basis 

(the Bell measurement) (I'L^)} = (8) and = 

{I,Z,X,ZX} is a set of operations to be applied conditional on the 
measurement outcome specified by k. Note that in case of n ^ Zd= 1, we have 
to repeat the teleportation protocol to transmit a quantum state between the 
nodes via the neighboring nodes. Thus all the Bell pairs corresponding to the 
vertical edges between l and n are consumed for performing teleportation. 


Therefore, Ci. nr .n is successfully applied on the control 
qubits at nodes vij and v rn j and the target qubit at node 
v n j by FOCC assisted by the Bell pairs corresponding to the 
vertical edges Sj between nodes Vy,j and v rn] . 


B. A conversion protocol 

We present a protocol to convert a given (k, A'(-cluster 
network into quantum circuits. First (step 1 to step 3), we 
construct quantum circuits of unitary operations that are im- 
plementable on qubits in a set of vertically aligned nodes VJ 
by FOCC assisted by the Bell pairs given for a set of vertically 
aligned edges Sj for a certain j. Then (step 4), we repeat the 
procedure given by the first part (step 1 to step 3) for different 
j of 1 < j < N. 

The conversion protocol is specified as follows: 
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1) Draw k horizontal wire segments where each of the wire 
segments corresponds to a set of qubits at vertically 
aligned nodes V' . 

2) Symbols representing two-qubit controlled unitary oper¬ 
ations Ci- n or three-qubit fully controlled unitary opera¬ 
tions Cz, m;rl are added on the horizontal wire segments 
according to the following rules. 

• To represent C; ;ra , draw a black dot representing 
a control qubit on the Yth wire, draw a vertical 
segment from the dot to the n-th wire segment and 
draw a box representing a target unitary operation 
on the n-th wire segment at the end of the vertical 
segment. Write index l at the side of the vertical 
segment between the horizontal wire segments. An 
example is shown in Fig. [6] i). 

• To represent Ci tm - n , draw two black dots repre¬ 
senting control qubits on the /-th and m-th wire 
segments, draw vertical segments from each dot 
to the n-th wire and draw a box representing an 
arbitrary target unitary operation on the 7i-th wire 
segment at the end of the vertical segment. Write 
indices l and m at the sides of the vertical segment 
between the horizontal wire segments. An example 
is shown in Fig. [6] ii) 

• Multiple gates of Ci ;n or C^ m . n can be added as 
long as there are only one type of index appearing 
between the horizontal wire segments and no target 
unitary operation represented by a box is inserted 
between two black dots on a horizontal wire seg¬ 
ment. An example of possible circuits generated in 
this protocol is shown in Fig. [6] iii). We also give 
an example of circuits that do not follow the rule in 
Fig. [6] iv). 

3) Arbitrary single qubit unitary operations represented by 
boxes are inserted between before and after the sequence 
of Ci ]n and Ctm-.n (but not during the sequence) ob¬ 
tained by step 2. 

4) Repeat step 1 to step 3 for each 1 < j < N and connect 
all the i-th horizontal wire segments. 

In Appendix A, we show that a unitary operation repre¬ 
sented by the quantum circuit obtained by step 1 to step 
3 of the conversion protocol is implementable in a set of 
vertically aligned nodes VJ, namely, it is implementable by 
LOCC assisted by (k — 1) Bell pairs corresponding to a set 
of vertically aligned edges Sj. As examples, quantum circuits 
converted from the (2, 3)-cluster and (3, 2)-cluster networks 
are shown in Fig. [7] 

Our conversion method generates infinitely many quantum 
circuits in general. However for special cluster networks, stan¬ 
dard forms of quantum circuits can be obtained. In Appendix 
B, we show that any converted circuit obtained from a (2,3)- 
cluster network can be simulated by the circuit presented in 
Fig. m i), and any converted circuit obtained from a (3, 2)- 
cluster network can be simulated by the circuit presented in 

FigHHi). 



Fig. 6 . i) A symbol representing Cz-i. ii) A symbol representing C\^%. iii) 
An example of circuits generated in step 2 of the conversion protocol. The 
index in the upper region is 1 , that of index in the middle region is 1 and 
that of index in the lower region is 4. iv) This conversion is forbidden since 
there is a target unitary operation inserted between two black dots representing 
controlled qubits. 


ii) 


Vi 

Cl, 3 


C \-2 




Vi 

Cl , 3 


C 2-1 


Vs 


i 1 i 1 i 

i 3 i 3 ! 


Vi 
GY, GY 


V 2 


HTh-p 

-m-Q 


Vi 


C2-I 


{pep—q> 

Hi—dm—ckh- 


V 2 


Fig. 7. Arbitrary single qubit unitary operations are represented by boxes, i) 
An example of converted quantum circuits from the (2, 3)-cluster network. 
It is obtained by connecting three segments of circuits generated in step 
1 to step 3 of the protocol corresponding to VJ', V 2 and V 3 . It consists 
of two-qubit controlled unitary operations defined by C7 Z;n = |0)(0| z 0 
+ | 0 )( 0 | z 0 Unl where l denotes the wire segment of the control 
qubit and u$ are arbitrary single qubit unitary operations on the n-th 
qubit, ii) An example of converted quantum circuits from the (3, 2)-cluster 
network. It consists of three-qubit fully controlled unitary operations defined 

by Q,m;n = |00>(00| Z , m 0 + |01>(01| Z , m 0 + |10)(10| Z , m 0 

^n 10 ) + |ll)(ll|z m 0 Un^\ where l,m denotes the wire segments of the 
’ c 

control qubits and are arbitrary single qubit unitary operations on the n- 
th qubit, iii) Another example of converted quantum circuits obtainable from 
the (3, 2)-cluster network. 


V. Implementability of unitary operations over 

THE BUTTERFLY AND GRAIL NETWORKS 

For classical network coding, it has been shown that there 
exists a network coding protocol over a 2-pair network, which 
has two input nodes and two output nodes, if and only if 
the network has at least one of the butterfly, grail or identity 
induced substructures ED, EH. Thus any classical network 
coding protocol over a 2-pair network can be reduced into 
a combination of protocols over the butterfly, grail or iden¬ 
tity networks, and these networks are fundamental primitive 




























































Fig. 8 . The nodes ij, i 2 > 01, 02, n 1 and n 2 of the butterfly network corre¬ 
spond to the nodes Di, 1, 1)3,1, 1)1,2, 1)3,2, 1)2,1 and 1)2,2 of a (3, 2 )-cluster 
network each other. The two-qubit unitary operation U g i 0 t, a i(x,y,z ) = 
e t(xX®x+yY®Y+zZ<giZ) ; s j m pl e mentable over a (3, 2 )-cluster network 
by fixing an input state at the node 1)2,1 at 10 ), performing an appropriate a 
three-qubit unitary operation U 3 and performing an appropriate LOCC map 
r consisting of a measurement on the qubit at the output node 1)2,2 and the 
conditional operations on the other output nodes 1)1,2 and 1)3,2 depending on 
the measurement outcome. 

networks for classical network coding. As a first step to 
investigate the implementability of quantum computation over 
general 2 -pair quantum networks, we investigate the imple¬ 
mentability of two-qubit unitary operations over the butterfly 
and grail networks in this section by using the method for 
converting a (fc, TV)-cluster network into quantum networks 
introduced in the previous section. 

We consider a two-qubit unitary operation U given in the 
form of the Kraus-Cirac decomposition represented by Eq. ©. 
Since it is trivial that the single-qubit unitary operations v and 
v' are implementable at the input nodes and w and w' are 
implementable at the output nodes, we just need to analyze 
the implementability of the two-qubit global unitary part 
U global (%, V , z) given by Eq. © over the butterfly and grail 
networks. Owing to the implementability of a three-qubit fully 
controlled unitary operation over the butterfly network, we 
have discovered a protocol for implementing U g i 0 b a i{x,y, z) 
for arbitrary x, y, z as presented in the constructive proof of 
the following Theorem. 

Theorem 1. Any two-qubit unitary operation is deterministi¬ 
cally implementable over the butterfly network. 

Proof. For the implementability of U g i 0 bai{x, y, z) over the 
butterfly network represented by the left hand side of 
Fig. 1 we consider a (3,2)-cluster network represented 
by the right hand side of Fig. [ 8 ] by assigning the nodes 
{*i,nr,* 2 ,Oi,ri 2 , 02 } of the butterfly network to the nodes 
1,^3,1,^1,2,^2,2,^3,2} Of the (3,2)-cluster network, 
respectively. In this assignment, the correspondence of the 
edges of the butterfly network and the horizontal and vertical 
sets of edges /Ci,<Si,«S 2 of the (3, 2 )-cluster network is given 
by 

{Ei, E 5 , E 3 } «-»■ /Ci, 

{E 2 ,E^} o <Si, 

{E 6 ,E 7 }, o 5a. (17) 

Thus any two-qubit unitary operation is deterministically im¬ 
plementable over the butterfly network if any U g i 0 bai{x , y , z) 


in the form of Eq. © is deterministically implementable over 
the (3, 2)-cluster network where input states are given at nodes 
V 11 and 1 ) 3,1 and output states are obtained at nodes iq 2 and 
v 3 , 2 , since the topology of the butterfly network is the same 
as that of the (3,2)-cluster network. 

We construct a protocol implementing two-qubit unitary 
U g iobai{x,y,z) by setting a fixed input state at node 1 ) 2,1 and 
arbitrary two-qubit input state at nodes 1 ) 1,1 and 1 ) 3,1 as a 
three-qubit input state at input nodes I = { 1 ) 1 , 1 , 1 ) 2 , 1 , 1 ) 3 , 1 }, 
and implementing a three-qubit unitary operation denoted by 
U 3 over the (3,2)-cluster network followed by an FOCC map 
denoted by T performed at output nodes O = { 1 ) 1 , 2 , 1 ) 2 , 2 , 1 ) 3.2 }■ 
Recall that a unitary operation represented by the quantum 
circuit shown in Fig. [7] ii) is implementable over the (3,2)- 
cluster network. That is, two three-qubit fully controlled 
unitary operations Ci, 3 ; 2 are implementable, one at nodes I 
and another at nodes O. The following protocol shows that 
by choosing appropriate parameters for one of the three-qubit 
fully controlled unitary operations and one of single-qubit lo¬ 
cal unitary operations in U 3 , we can implement U g i 0 b a i ( x , y, z) 
with arbitrary x,y,z. 

The protocol for implementing U g i 0 bai{x,y, z): 

1) An arbitrary two-qubit input state p is given for qubits 
at input nodes iq,i and 1 ) 3,1 and a fixed input state |0) 
is prepared for the qubit at node 1 ) 2 , 1 . 

2) Implement t / 3 of which the quantum circuit representa¬ 
tion is given by the left shaded part of Fig. [9] over the 
(3,2)-cluster network. 

a) All single-qubit unitary operations appearing in the 
circuit representation of U 3 are trivially performed 
at each node. 

b) The first fully controlled unitary operation im¬ 
plemented at input nodes I using the Bell pairs 
represented by vertical edges 5) is given by 

Ci,3;2({t4 ab) }a, 6 = 0 , 1 ) where i4 00) = ui n) = I 
and it! 01 ) = t4! 0) = Z. 

c) To transmit a qubit state from input nodes Dj,i 
to output node Dj , 2 for i = 1,2,3, quantum 
teleportation is performed for each i by using the 
Bell pair represented by a horizontal edge in /Ci. 

d) The second fully controlled unitary operation im¬ 
plemented at output nodes O contains parameters y 
and z and is given by C[ 3 . 2 ({w^} a ,b=o,i) where 

id! 00 ) = «4 U) = e i(z_!/ )|0)(0| - ie^ z+v) |1)(1|, 
id! 01 ) = «4 10) = e -i C-D)|o)(o| - ie _i ( 2+y )|l)(l|. 

e) After implementing C! 3 ; 2 ({ii)! Qfc) } a) b = o,i), a 
single-qubit unitary operation parameterized by x 
given by 

U ^ = Tf ( e“ ie~ ix ) ( ' 18) 

is performed at node d 2 ,2 £ O. 

3) Perform an FOCC map T at output nodes O of which 
the quantum circuit representation is given by the right 
shaded part of Fig. [9] The map T consists of the 
following three steps. 
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Uz r 



lo> 


Fig. 9. A quantum circuit representation of a three-qubit unitary operation 
U 3 (the left shaded part) and an LOCC map T (the right shaded part) used in 
a protocol for implementing a two-qubit unitary operation U g i 0 b a i(x , y , z) = 
e %(xX®x+yY®Y+zZ®z) on t fo e fi rst anc j qubits. The input state of the 
second qubit is fixed in |0). The single-qubit unitary operations represented 
by boxes labelled by H and X are given by H = (|0)(0| + |0)(11 +|1)(0| — 
|l)<l|)/x/2, u(x) = H(e ix |0)(0| - ie-^|l)<l|) and X = |0)(1| + |1)(0|, 
respectively. The target single-qubit unitary operations of the first three- 
qubit fully controlled unitary operation Ci^- 2 ({u^ b ^} a b= 01 ) arc given 
by Un °^ = 1 and = Z. The target single-qubit 

unitary operations of the second three-qubit fully controlled unitary operation 
C ’l,3;2({"'^ l6> }a,i=0,l) && g iven b Y w n^ = 10) (0| - 

ie»U’+w)|i)(i| and «4 01) = «4 10) = |0>(0| - ie~ i( - z+v ') |1)<1|. 

The half circle symbol represents a projective measurement in the com¬ 
putational basis {|&)(&|}fc = o i- The single-qubit unitary operations (boxes) 
connected to the measurement symbol by dotted lines represent conditional 
unitary operations performed only if the measurement result is k = 1 and do 
nothing (or perform I) if k = 0. 


a) Perform a projective measurement on the 
qubit at node V 2 2 in the computational basis 
{|0)(0|,|1)(1|}. 

b) Classically communicate the measurement out¬ 
come k £ { 0 , 1 } from node ?; 2 , 2 to i>i 2 and also 

to l> 3 j2 - 

c) If k = 1, perform a conditional operation X on 
output qubits at nodes r>i , 2 and 113 , 2 , otherwise do 
nothing. 

This protocol maps any input state p given at input nodes 

i>i,i and 113.-1 to 

U global (x,y,z) P Ul lobal (x,y,z) = r(U 3 (p ® |0)<0|)C/+) (19) 

at output nodes v\p and 113 2 where | 0 ) represents the fixed 
input state at node ti 2 ,i. See Appendix D for details of 
calculations. It is straightforward to translate the protocol 
over the (3, 2)-cluster network to a protocol to implement 
U g iobai(x , y, z) over the butterfly network by using the corre¬ 
spondence of vertices and edges. Thus, U g i 0 b a i{x,y : z) is de¬ 
terministically implementable over the butterfly network. □ 

In the implementation protocol of U g i 0 b a i(x, y, z) presented 
in the above proof, the first-three qubit fully controlled op¬ 
eration Ci, 3 ; 2 {{u^}a,b=o,i) where ul° 0) = = I and 

Un 1 ) = Un {Y) = Z can be decomposed into a sequence 
of two controlled-Z operations Ci- 2 {{u^ = I, u'.p = Z }) 
and € 3-2 = ({'tt 2 °' > = I; U'a 1 = %})■ This sequence of two- 
controlled Z operation can be implementable by consuming 
two Bell pairs corresponding to the edges E2 and £4. On the 
other hand, the second three-qubit fully controlled operation 
C[ t 3 ; 2 ({wi a6) }o, 6=0,1) where «ii° 0) = w £ n) = e l(z_y) | 0 )( 0 |- 
iJ( z +y)\l){l\ and w { u 1] = w { n 0) = e^ 2 "^ |0)(0| - 

ie~ t< ' z+y ' 1 11 )(1| cannot be decomposed into two two-qubit 



Fig. 10. The nodes m, n 2 , 01 , i 2 , 713 and 714 of the grail network correspond 
to the nodes «i 1 , r>i, 2 , Vi .3, V2 1, V2 .2 and V2 ..3 of a (2, 3)-cluster network, 
respectively. The set of all unitary operations implementable over the (2, Si- 
cluster network is also implementable over the grail network, since we can use 
the edges E\ and i? 2 for just teleporting qubits and the rest of the network 
forms the (2, 3)-cluster network, with which any two-qubit unitary operation 
is implementable. 


controlled unitary operations which are implementable by 
using a Bell pair for each. This is the reason why direct 
implementability of a three-qubit fully controlled unitary oper¬ 
ation by just consuming two vertical Bell pairs corresponding 
to Eq and £7 is the key for proving implementability of 
U g iobai(x,y,z) over the butterfly network. 

For the implementability of U g i 0 b a i(.x : y, z) over the grail 
network, we consider a ( 2 ,3)-cluster network by assigning 
the nodes {m,n 2 , 01, i 2 , 77 . 3 , 714 } of the grail network to the 
nodes {i>i,i, ni, 2 , 111 , 3 , ti 2 ,i, 112,2, 12,3} of the (2, 3)-cluster net¬ 
work, respectively (Fig. [IS. The (2,3)-cluster network can 
be converted to a quantum circuit containing three controlled- 
NOT operations and arbitrary single unitary operations that are 
inserted between the controlled-NOT operations. It is shown 
that any two-qubit unitary operations U g i 0 b a i{x,y, z) can be 
decomposed by three controlled-NOT gates and single unitary 
operations inserted between the controlled-NOT operations 
El. Thus any two-qubit unitary operation is deterministically 
implementable over the grail network. 

VI. The set of all implementable unitary 
OPERATIONS FOR k = 2,3 

In this section, we derive the condition for fc-qubit unitary 
operations to be implementable over a given cluster network. 
We show that our conversion method presented in Section III 
gives all implementable unitary operations over the (fc,7V)- 
cluster network for k = 2,3. 

Theorem 2. If i) a k-qubit unitary operation U is deter¬ 
ministically implementable over the (fc, N)-cluster network 
(k > 2, N > 1), then ii) the matrix representation of U in 
terms of the computational basis U M can be decomposed into 

U M =V 1 M V 2 M (20) 

where each V\ M is a 2 k by 2 k unitary matrix such that 

v,i 4 = ee- E ^ } ®4r 2) ®4r 3) 

ai— 0 a,2— 0 afc-i — 0 


(21) 
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where E < : n l '" n> and E^ are 2 by 2 complex matrices. 

To prove Theorem 2, we first prove Lemma [j] about a class 
of bipartite separable maps that preserves entanglement. A 
bipartite separable map T sep is a completely positive and trace 
preserving (CPTP) map whose Kraus operators are product as 
follows: 


r sep(PEF) = Y,( Ek 0 Fk)pEF(Ek 0 F k )\ (22) 

k 

where J2 k ( E k 0 E k)H E k 0 F k ) = I E 0 If- Since quantum 
network coding is equivalent to perform LOCC assisted by 
the resource state in our setting, we have to analyze multi¬ 
partite LOCC. However, the analysis of multipartite LOCC 
is extremely difficult. Thus, we analyze multipartite separable 
maps, which are much easier to analyze than LOCC due to 
their simple structure. Note that a set of separable maps is 
exactly larger than that of LOCC l27l . 

Let |V^ n ) = \ A i)\ B i) and I ^out) = 

75 XEo l a i)l & i)> where {|A) G Ha} and {| Bi) G U B } 
are orthonormal sets and {|oj) G H a } and {|h,) G Hb} are 
orthonormal bases. Note that dim('H 0 ) = dim(’Hb) = d but 
the dimensions of Ha and Hb can be higher than d, therefore 
{|Ai)} and do not need to form bases. 

Lemma 1. Let {E k G L (Ha '■ H a )},{F k G L (Hb ■ H b )} be 
sets of linear operators. If {E k 0 F k } satisfies 


Y J E l E k®FlF k = Iab (23) 

k 


and for all k. 


By using Eq. ( l26t and Eq. ( |271 >. we obtain 
tr (ll E k ]E k ® E k ] F^j 
= tr( E E ^ E ™® E ^ E k M ) +c 

\fce{fc|j>fc#o} / 

E P ktT ( E ^ E k ® (Ejf'Ej?)-') +e = d 2 

k£{k\p k ^0} 

E P* tT { E ^ E k ® (E^Ek)- 1 ) = d 2 -e, 

fce{febfe#o} 

(28) 

where e = tr (£ fce{fc | Pfe=0} E k'* E k 0 F k* F k') > 

We let P k = E^ Ejf be a d by d positive matrix and 
{A^, > 0|i = 0,1, • • • , d — 1} be the set of eigenvalues 
of P k . Then the eigenvalues of = Pjf 1 are 

= 0,1, • • • , d — 1} and the condition Eq. (l28l) is given 

E PkYKYjj = d 2 -e. (29) 

ke{k\p k ^0} i=0 j=0 A k 


Using the Cauchy-Schwarz inequality, we obtain 


d— 1 d— 1 -j 

= 

4=0 j =0 A k 


> 



(30) 


The equality holds if and only if Aj. = a 2 > 0 for all i. 
By using Eqs. d29t-d30t and the fact that {p k \pk f 0} is a 
probability distribution, we obtain for all k G {k\p k ^ 0}, 


E k 0 Pfc|4/j n ) — \JPk T out } (24) 


is satisfied, then for all {k\p k 0 }, 


3a > 0; P k = Ef f Pf = a 2 I d , (31) 

e = 0. (32) 

□ 


3a fc > 0, 3C/f G U(C d ), Pf = a k U™,F™ = , 

Oik 

(25) 

where Ej}' 1 and F^ 1 are d by d matrices such that (Ejfi)ij = 
(ai\Ek\Aj), {F^) it j = (bj\ F k \Bj), U(C d ) is the set of d by 
d unitary matrices and U M is the complex conjugate of U M . 


Proof. By straightforward calculation, we obtain 

Vfc, Pf(Pf) T = y/pih 
=> Vfc G {k\p k ± 0}, P fe M = y/pf((E^)~ 1 ) T , (26) 

and 

E(£f = I d 2. (27) 


Proof of Theorem 2. Denote by Hi q = ® k -iHi i and 
Ho q = ®i—\HOi the Hilbert spaces of k input qubits and 
k output qubits, respectively. By introducing another ancillary 
Hilbert space Hp at the input nodes t', j, denote the Hilbert 
space of k qubits by Hi> = ® k =l Hp. A joint state of k 
copies of the Bell pairs in Hi q 0 Hi' Q is denoted by 

i 1 ) : =-4 ei*>^h=0*=iI* + >^' 

v i 

where D = dim (Hi q ) = 2 k . If U G \J(Hi q : Hq q ) is 
deterministically implementable over a (k. iV)-cluster network 
for k > 2 and N > 1, it is possible to apply U on |I) and to 
transmit the resulting state to the output nodes. That is, there 
exists a LOCC map T such that 

^E r (l i )O1x Q ®l^}(^k)0|*)O1x i =|C/)(f/|, (33) 

h3 
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where |<f>) 7 j is the resource state of the (k, 3V)-cluster network 
and |Z7) is defined by 


Let 


\U) := ( U ® I)|I) G 7-Loq ® T~Li' Q - 


(34) 


By defining a map represented by the left hand side of 
Eq. m as r(|$)<$| w ) := £ J2i tj r(|*)(i|x 0 ® |*><*| W ) ® 
\i)(j\x' , where T' is also a LOCC map if we assume two 
qubits belonging to 77 /, and 77 r are in the same input node 
for all i. Since any LOCC maps are separable maps, there 
exists a separable map V sep satisfying 


Ke P m($\n) = \U)(U\, 



= U S i , j 

3 77 ft', 

,3 

(39) 

r^ 


(* 

= D 

(40a) 



(2 

< i < k - 

- 1) (40b) 



(* 

= k) 

(40c) 

^ K h 


0 

= D 

(41a) 

U *L 

® 77 r /2 

(2 

VI 

VI 

- 1) (41b) 

% KIn 


0 

= iV) 

(41c) 

m = ® 

fc /im 

2 = 1 ^, l> 

-Pm = 

®J =1 ® jL 

•2 -A^ and 




and 


= E^ka 0 ^’™ 

ai —0 
1 1 

a™ = E 

a± —0 a 2 —0 

(2 < i < k- 1 ) 

ai =0 


(ai,a 2 ),^ 


(43) 

(44) 

(45) 


(35) 


where 77 


ft 0 ’™ G L(^ ni : 77 /( ), G L(W ni : 

77/') and E ( k a f’ m G 1,(77*^ : 77/'). Thus, can be 
decomposed into 


(oi,a 2 ),m 


if 17 is deterministically implementable over a (7, A')-duster 
network. Since r' ep is a map from a pure state to a pure state, 
the action of r' represented by Ea. (l35l) can be equivalently 
given by the existence of a set of linear operators (the Kraus 
operators) {.4™-}m for each node Vij and a probability 
distribution {p m } such that 

Vm; ®? = i ®jLi = VP^\U), (36) 

&) = I, (37) 

m 

where 

G L(H u<il : 77/') 

(1 < i < k), 

A™ G L (U ViJ : C) 

(1 < i < k,2 < j < N - 1), 

G L (H Vi , N :H 0i ) 

(1 < i < k), (38) 

and / H Vi . is the Hilbert space of qubits of the resource state 
at node Vij defined by 


Vi 


M 
1 ,m 


E E Ci 

ai,-" i a k — l —0 




Note that we identify a linear operation and its matrix repre¬ 
sentation in the computational basis, e.g., ' t,m is a 2 by 2 
complex matrix. 

Next, letting E m = ®f =1 <g >? =1 A™-, F m = ®£ =1 <g )^ 3 
and using Lemma [T| we obtain for all m G {m|p m ^ 0}, 


3ct2,m 

> 

0,3L 2 M m GU(C fl ); =« 2 

v M 

,m v 2 ,m 

where E™ 

is 

a 77 x 77 matrix such that 


(C)«, 6 

= 

kk(k=i ®? =1 ^)l^ 6 )s ; . 

i,S:, 2 ,if; 

and 




l^i)s 

7,i 

, S t,2’ K l,i’ K Z,2 = 0 i=l 1 |^ >+ )s , Ji 




=i 1 |^ + )s t 1 , 21 sf +1>2 ®i=t 





■^, 2 - 


applying Lemma Q] we obtain for all m G {?n|p m ^ 0}, 

3a hm > 0, AV^ G U(C D ); E™ = (42) 

where U(C D ) is the set of D by D unitary matrices and E 
is a D by D matrix satisfying 


I— ®i=l k + )s, 1 1 ,S ? +11 ® L - 


Let 


^2 = E< a ik 2 0 ^L2 l),m (48) 

CL\ —0 

= E X>k>k a 0 ^ ),m 

ai —0 a,2=0 

(2 < i < k - 1) (49) 

1 

^2 = E^l^, 2 0 4?2 l),m > (50) 

ai —0 

where E[f’ m G L(H K i 2 ® 77 ^ 2 : C), G 

L(W x i 2 <3 77/^ : C) and E$’ m GL^®^ : C). 
By straightforward calculation, 17/^ are shown to be decom¬ 
posed into 

V 2 M m = V ft, £ E^ m ®E’ 2 ^ m ®... 

a- 1 ,-" ,afc_i =0 

®^-i; 2 2 ’ a '" l) ’ ro 0 (51) 
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where E^ 1 ^’™, anc j E^ff' 1 '"’ are 2 x 2 complex 

matrices. 

Iterating this procedure, we obtain for all m £ {m\p m ^ 0}, 


3a > 0, 3W m £ U(C D ); E% = aW M ,F™ = 

where W M and W M can be decomposed into 



(52) 


W M 


= V, m V 2 m 
W M = 


...V M 

V N -1 


(53) 

(54) 


and V\ M 


V 1 

,afc_i=0 1 ,i 


e: 


(ai,a 2 ) 


E k-Ci' ak ~ l) ® E k7~ l] e U(C D ). U M can be decomposed 
into the form of La. (Eoli since V t M and can be decom¬ 
posed into the form of Eg. (EH. 


□ 

In the case of the (2,7V)-cluster networks, which we call 
iV-bridge ladder networks, Vi is locally unitarily equivalent 
to the two-qubit controlled unitary operation since its operator 
Schmidt rank is 2 lE6l . Thus, statements i) and ii) of Theorem 
2 are equivalent since a sequence of N two-qubit controlled 
unitary operations is implementable by the converted circuit 
presented in Fig. [7] Then we obtain the following theorem for 
the ladder networks. 

Theorem 3. A unitary operation U is deterministically im¬ 
plementable over the N-bridge ladder network if and only if 
KC#(17) < N. 

This theorem is proven by using the following lemma relat¬ 
ing the Kraus-Cirac number of a two-qubit unitary operation 
and the decomposition of the unitary operation into controlled 
unitary operations shown in Ill- 

Lemma 2. Consider a set of two-qubit unitary operations U c 
that is locally unitarily equivalent to a controlled unitary op¬ 
eration. The decomposition of a unitary operation U £ SU(4) 
into a shortest sequence of two-qubit unitary operations in U c 
depends on the Kraus-Cirac number KC#(U ) of U as 

{U £ S , (7(4)|KC#(/7) < 1} = {U\U £ U c } 

{U £ 577(4) |KC#((7) < 2 } = {UV\U,VeU c } 

{U £ 577(4) |KC#(7/) < 3} = {UVW\U, V,W € U c j. 


VII. Probabilistic implementation of unitary 

OPERATIONS 

It is interesting to know whether there exists a task that 
is not achievable by classical network coding but the cor¬ 
responding task is achievable in a quantum setting. We can 
give a negative result in the case of a ( 2 , 2 )-cluster network. 
There is no classical network coding protocol to send single 
bits from v\p to 1 * 2,2 and from 1 * 2,1 to 1 * 1,2 over a ( 2 , 2 )- 
cluster network since there is no butterfly, grail or identity 
substructure. This task corresponds to implementing a SWAP 
operation in quantum network coding. By using Theorem 
3, we see that a SWAP operation is not deterministically 
implementable over a ( 2 , 2 )-cluster network, which is a 2 - 
bridge ladder network, since the Kraus-Cirac number of the 
SWAP operation is 3. 

Unitary operations are deterministic maps by definition, but 
we consider the less restricted situation where the action of 
the unitary operations are implemented only when we can 
post-select the preferable probabilistic event. This corresponds 
to requiring the implementing the action of a unitary opera¬ 
tion only when certain measurement outcomes in a LOCC 
protocol are probabilistically obtained. A formal definition 
of the probabilistic implementation of a unitary operation 
is given by Definition 3 by changing LOCC to stochastic 
LOCC (SLOCC). In this section, we first characterize all 
the unitary operations that are probabilistically implementable 
over cluster networks. Then, we show that a SWAP operation 
is not implementable even probabilistically. 

Theorem 4. A k-qubit unitary operation U is probabilistically 
implementable over the {k, N)-cluster network (k > 2,N > 
1) if and only if the matrix representation ofU in terms of the 
computational basis U M can be decomposed into 

U M = pM p M ... F M } ( 55 ) 

where each F t M is a 2 k by 2 k complex matrix that can be 
decomposed in the same way as Eq. ( 1271 ) 

Proof. Similar to the case of deterministic implementation, 
we consider applying U £ U (TLx Q '■ T~Lo Q ) on a part 
of k maximally entangled states |I) £ TLx Q <8> Then 

U is probabilistically implementable over the ( k , A'(-cluster 
network (k > 2, N > 1) if and only if there exist a stochastic 
LOCC (SLOCC) map T" and non-zero probability p > 0 such 
that 

r"(|$X$k) =p\u){u\, ( 56 ) 


Proof of Theorem 3. Since KC#(U) is less than or equal 
to N if and only if U can be decomposed into N two-qubit 
controlled unitary operations as shown in Lemma [2] and N 
two-qubit controlled unitary operations are deterministically 
implementable over A’-bridge ladder network, Theorem 3 is 
straightforwardly shown. 

□ 

We also show that statements i) and ii) of Theorem 2 
are equivalent in the case of the (3, iV)-cluster networks in 
Appendix C. 


where | < I>)-k is the resource state of the (k, W)-cluster network 
and | U) £ Tio Q <S> T~l-x' Q is defined by Eq. (E4t . Eq. ( f56t 
is equivalent to the statement that there exist a set of linear 
operators {A l] \ and non-zero probability p > 0 such that 

= (57) 

The conditions of {Ajj} given by Eq. (l57l i is similar to 
the conditions of Kraus operators {A™ } m given by Eq. ( l36l > 
presented in the proof of Theorem 2. The index m is dropped 
in Eq. (l57l > since the map we consider is SLOCC instead of 
LOCC considered in Theorem 2. By taking the correspondence 
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between A, j and ,4”), we obtain a decomposition of the form 
presented in Eq. ([55l i. 

□ 

Lemma 3. A SWAP operation U swap := |00)(00| + |01)(10| + 
110 ) (011 + 111 ) (111 is not probabilistically implementable over 
the 2-bridge ladder network. 


Proof. By using Theorem 4, the SWAP operation is probabilis¬ 
tically implementable over the ( 2 , 2 )-cluster network ( 2 -bridge 
ladder network) if and only if there exist linear operations 


P, Q S L(77i <g> 77 2 ) and Ef) € L (Hf satisfying 


Uswap — UQi 

(58) 

P = 

(59) 

Q = Ei°lQE$+EW®EW t 

(60) 

where 77,; = C 2 . Since P and Q can be decomposed into 

Ea.(l59l> and Eq.(f60t. we can derive 


Op#i(P) < 2, 

(61) 

Op#i(Q) < 2. 

(62) 

Since Op #l(U swap ) = 4, Op#?(P) = Op #\(Q) 

= 2. In 


[|28l . it is shown that if Op#^(P) = 2 and P is invertible, 
Op#x(P _ 1 ) = 2. Thus, the SWAP operation is probabilisti¬ 
cally implementable if and only if there exist linear operations 


P, Q S L(77i ® P 2 ) satisfying 

Q = U swap P , (63) 

Op#i(P) = 2, rank(P) = 4 (64) 

Op#i(Q) = 2, rank(Q) = 4. (65) 

In general, we can regard P as a matrix representation of a 
four qubit pure state | 3 >)i, 2 , 3 , 4 ; 

4 

P = 5Z(* |3,4|$)l,2,3,4(*|l,2- (66) 

*=1 

Then, the following correspondences are obtained, 

rank(P) = 4 SCH#^^)) = 4, (67) 

Op #l(P) = 2 Sch#^3(|$)) = 2, (68) 

Op#?(P su;op P) = 2 o Sch#^(|$)) = 2, (69) 


where SCHtt^d^)) is a Schmidt number in terms of a 
partition between qubit 1, 2 and qubit 3,4. We show that there 
is no four qubit state simultaneously satisfying Eqs. 
and (l69l > in Appendix E. □ 

We can apply Theorem 2 and 4 to a slightly extended cluster 
network, a cluster network with loops. We show the definition 
in Appendix F. 


circuit representations of unitary operations implementable 
over a given cluster network. For the (k, A')-cluster networks 
of k = 2,3, we have shown that our method provides all 
implementable unitary operations over the cluster network. 
The proof is based on the existence of the standard form of 
the converted quantum circuit and the equivalence of a set 
of unitary operations represented by the standard form and 
decomposed into the form given by Eq. <ED. The proof also 
suggests that statements i) and ii) of Theorem 2 are equivalent 
for k = 2,3. For k > 4, whether our method provides 
all implementable unitary operations or not is still an open 
problem since the standard form is not known. 

As a first step to finding the fundamental primitive networks 
of network coding for quantum settings, we have shown that 
both of the butterfly and grail networks are sufficient resources 
for implementing arbitrary two-qubit unitary operations, mean¬ 
while the ( 2 , 2 )-cluster network is not sufficient to implement 
arbitrary two-qubit unitary operations even probabilistically. 
To prove this, we have shown necessary and sufficient con¬ 
ditions of probabilistically implementable unitary operations 
presented in Theorem 4. There are two differences in Theorem 
2 and Theorem 4. First, we have shown that U can be 
decomposed into a particular form represented by Eq. <[20b 
ifU is deterministically implementable in Theorem 2 and that 
U can be decomposed into a particular form represented by 
Eq. (1551) if and only if U is probabilistically implementable 
in Theorem 4. Second, each factor in Eq. (l55l ) can 
be a non-unitary complex matrix while each factor V) M in 
Eq. (l20l) must be a unitary matrix. The existence of unitary 
operations only probabilistically implementable (with less than 
unit probability) is also left as an open question. 
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VIII. Concluding remarks 

We have investigated the implementability of fc-qubit uni¬ 
tary operations over (k, A - (-cluster networks to apply the idea 
of network coding for distributed quantum computation where 
the inputs and outputs of quantum computation are given in all 
separated nodes and quantum communication between nodes 
is restricted. We have presented a method to obtain quantum 
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Fig. 11. An example of a converted quantum circuit obtained by step 2 of 
the conversion protocol. 


Appendix A 

LOCC IMPLEMENTATION OF CONVERTED QUANTUM 
CIRCUITS 

We have shown a protocol to implement a three-qubit fully 
controlled unitary operation in a set of vertically aligned 
nodes VJ. In some cases, we can implement more than 
one three-qubit or two-qubit controlled unitary operations in 
parallel using the same resource. We show how a sequence of 
controlled unitary operations represented by converted circuits 
can be implemented by LOCC assisted by the resource state 
given by a collection of (k — 1) Bell pairs corresponding to a 
set of vertically aligned edges Sj in this appendix. 

We introduce a new notation for controlled unitary opera¬ 
tions for simplifying and unifying descriptions of two-qubit 
and three-qubit controlled unitary operations. We represent a 
two-qubit controlled unitary operation that is controlled by the 
7-th qubit and targets the j-th qubit as 

(70) 

and a three-qubit fully controlled unitary operation that is 
controlled by the z-th and j-th qubit and targets the /.:-th qubit 

as 

(i,r,k). (71) 

Note that we represented () as Ci-j and ( i,j;k ) as 
Cjjjfc in the previous sections. Let G = {g n } be a sequence 
of controlled unitary operations that is added in step 2 of 
the conversion protocol. For example, the converted circuit 
represented by Fig. Qj]is described by a sequence 


51 

= ( 1 , 1 ; 2 ) 

(72) 

52 

= (4,4; 2) 

(73) 

53 

= (1,4; 2) 

(74) 

54 

= (4,4; 5) 

(75) 

55 

= (4,4; 3) 

(76) 

56 

= (5, 5; 6 ) 

(77) 

57 

= (4,4; 5). 

(78) 


Let Ci be a set of controlled unitary operations that is 
controlled the 7-th qubit: 


(79) 


For example, for G defined by Eqs. (1721)- (T78l) . 

Ci = { 91 , 93 } (80) 

C4 = {32,53,34,35,37} (81) 

C 5 = {35} (82) 

C 2 = C 3 = Cq = 0. (83) 


Define the range of C, ^ 0 as 


range(Ci) = (min{? 

!, min{(a, 6; c) £ CJ}, 

C 


max {7 

, max{(a, 6; c) £ Ci}}). 

C 

(84) 

For example, for Ci defined 

by Eqs. S-®, 


range (Ci) 

= ( 1 , 2 ) 

(85) 

ranged) 

= (2,5) 

( 86 ) 

range (C 5 ) 

= (5,6). 

(87) 


All the controlled unitary operations in G are implementable 
by using the following protocol. 

The protocol for implementing a sequence of controlled 
unitary operation in G : 

1) For applying gates in Ci, we create an ancillary qubit 
state entangled to the v'-th qubit state by preparing an 
ancillary qubit in |0) and applying a controlled-NOT 
operation where the ancillary qubit is the target qubit 
of a controlled-NOT operation. Then the ancillary qubit 
state is sent from the z-th node Vij £ VJ to the target 
node by using teleportation. If several different target 
qubits are included in Ci, create another ancillary qubit 
by the same method at a target node, keep one of the 
ancillary qubits at the target node and send the other 
to the next target node. We consume ru Bell pairs to 
teleport ancillary qubit states to the target nodes, where 
Hi = b — a and range(Ci) = (a, b). Since there is no 
overlap between ranges of Ci and there is no target 
unitary operation inserted between control qubits, we 
can teleport all the ancillary qubit states entangled to the 
control states to all the target nodes by just consuming 
[k — 1) Bell pairs. 

2) We apply all the controlled unitary operations in G in 
the target nodes by using the teleported ancillary qubit 
states entangled to the control qubit states as the control 
qubits. 

3) We decouple the ancillary qubit states by performing 
the projective measurements on the ancillary qubits in 
the target nodes and apply correction unitary operations 
in the control nodes depending on the measurement 
outcomes. 


Appendix B 

Converted circuit of (2, TV) and (3, 7V)-cluster 

NETWORK 

First, we prove that any converted circuits of a (2,7V)- 
cluster network can be simulated by a circuit consisting of a 
sequence of TV two-qubit unitary operations and local unitary 
operations. In this case, only two-qubit unitary operations 
( 1 , 1 ; 2 ) or ( 2 , 2 ; 1 ) can be added in step 2 of the conversion 


Ci = {(a, b; c) £ G; a = i V b = i}. 
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i) ii) iii) 



Fig. 12. The six classes of converted quantum circuits obtained by step 2 of 
the conversion protocol of a (3, TV)-cluster network. 


protocol. Since applying the gate (1,1; 2) sequentially for 
k £ N times can be simulated by just one use of gate (1,1; 2) 
and gate (2, 2; 1) can be simulated by one use of gate (1,1; 2) 
and additional local unitary operations, any circuits generated 
in step 2 of the conversion protocol can be simulated by one 
use of (1,1; 2) and local unitary operations. 

Next, we prove that any converted circuits of a (3, N)- 
cluster network can be simulated by the circuit of a sequence 
of N three-qubit fully controlled unitary operations given in 
the form of 

|00)(00|i,3 0 U ( 2 00) + |01}(01|i,3 0 4 01) 

+|10)(10|i,3 0 U ( 2 10) + 111)<111 1,3 0 4 n) (88) 

and local unitary operations. In step 2 of the conversion 
protocol, every converted circuits can be simulated by six 
classes of circuits shown in Fig. 12. 

In the following, we show that all of these six classes (from 
class i) to class vi) represented in Fig. 12) can be simulated 
by a three-qubit fully controlled unitary operation and local 
unitary operations by investigating each class, 

i) A unitary operation obtained by circuit i) is given by 

|0)(0|i 0 u 2 0 U3°) + 11)(11 1 0 u 2 ^ 0 U 3 1 ' 

== |0)(0|i 0I 2 0 I 3 + 11)(11 1 0 0 u^u^ f (89) 

where is a one-qubit unitary operation and = rep¬ 
resents local unitary equivalence. Diagonalize 


3 0)t as 





M 2 1)u 2 0)t = 

= V2 1 

( e iSl 

l 0 

e i) 4 

(90) 

4 1) ' i 4 0)t = 

= v 3 1 

/ e ie 3 

l 0 

A «3- 

(91) 


Since the right-hand side of Eq. ( l89t is locally unitarily 
equivalent to a diagonal unitary operation in the compu¬ 
tational basis, this circuit can be simulated by a three- 
qubit fully controlled unitary operation and local unitary 
operations. 


ii) In circuit ii), the two-qubit controlled unitary operation 
(2, 2; 3) can be decomposed into 

|O)(O| 2 0 4° ) + |1)(1|2 0 4 1) 

= v 3 (jo)(o| 2 0 1 3 + |i)(i| 2 0 (^ e 0 A)) 44 0) 

= (I2 0 V 3 ) 

((0 A)® |0)(0|3+ (o A)® |i)(i|s ) 

(I 2 0 Ugttg 0 ^), (92) 

where v 3 is a unitary operation that diagonalizes 
?4 1) t4 0)t . Thus, this circuit is locally unitarily equivalent 
to a three-qubit fully controlled unitary operation. 

iii) Circuit iii) consists of just a three-qubit fully controlled 
unitary operation. 

iv) Circuit iv) can be simulated by a three-qubit fully con¬ 
trolled unitary operation and local unitary operations 
since we can diagonalize a unitary operation obtained by 
the circuit in the same way as circuit i). 

v) In the same way as circuit ii), circuit v) is locally uni¬ 
tarily equivalent to a three-qubit fully controlled unitary 
operation. 

vi) In the same way as circuit i), circuit vi) is locally uni¬ 
tarily equivalent to a three-qubit fully controlled unitary 
operation. 


Appendix C 

Two CONDITIONS IN THEOREM 2 ARE EQUIVALENT IN THE 
CASE OF THE ( 3 , iV)-CLUSTER NETWORKS 

For k = 3 , the 2 k by 2 fc unitary matrix V t M in Theorem 2 
is written by 


V M = E (°) 


l,i 

+E- 


+E\ 


+E- 


1 E. 


T ■ 


1 E. 


(o) 

' 3 ,i 


(0) 
1 ,i 


( 1 ) 


( 1 ) 
1 ,i 


E?f 


E$ 

4 ? 

e£1 


(93) 


By using the result on local unitary equivalence of unitary 
operations with operator Schmidt rank 2 obtained by Cohen 
and Yu l29l (Theorem 1 of fl29l ), we have 


V™ W |0)(0|a 0 W^c + |1)(1U 0 (94) 

= ^AB 0 |0)(0|c + W^B ® |1)(1|C) (95) 

where W^ c and W^ B are unitary matrices, = J represents a 
locally unitarily equivalence and we identify a three-qubit uni¬ 
tary operation on TLa 0 T~Lb ® 'He as i ts matrix representation 
V. { M . Thus, it is shown that 

V\ M L = |00)(00| J 4c 0 + |01) (01 |^4c 0 Wg lS> 

+|10)(10Uc ® W { B 10) + |ll)(llUc <8 wg 1} , 

(96) 

(ij) 

where W B is a 2 by 2 unitary matrix. Statements i) and ii) 
of Theorem 2 of the main text are equivalent in the case of 
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(3, iV)-cluster networks since l/ iW is a fully controlled three- 
qubit unitary operation and N fully controlled three-qubit 
unitary operations are implementable by a converted circuit 
of the (3, lV)-cluster networks. 


Appendix D 

A NETWORK CODING PROTOCOL FOR THE BUTTERFLY 
NETWORK IMPLEMENTING ARBITRARY TWO-QUBIT 
UNITARY OPERATIONS 

We show that the quantum circuit presented in Fig. [9] 
implements a two-qubit global unitary U g i 0 b a i{x,y, z) given 
by Eq.© for arbitrary parameters x,y,z £ K. U g i 0 b a i(x , y, z) 
can be decomposed into 

u global (x,y,z) = J2 X ^ U) )^ U) \ (97) 

3 

by using its eigenvalues {Aj}^ and eigenvectors {|\f , ( : 0)}j 
such that 


_ e i(x-y+z) 

,Ai = e i (- x +»+*), 

(98) 

— e i( x +y-z) 

', As = e i( - x - y ~ z \ 

(99) 

l«p( 0) ) = 

±(\00) + \U)), 

(100) 

|T (1) ) = 

±(\00)-\U)), 

(101) 

|«p(2)) = 

-L(\01) + \10)), 

(102) 

|q>( 3 )) = 

-^(|01)-|10)). 

(103) 


Thus, in order to show an arbitrary input state | \f/} is trans¬ 
formed into Uglobal (x, y, z) |\I/) through the quantum circuit, 
it is sufficient to show that the eigenvectors {|'k^))} J are 
transformed into {Aj | d/tj 1 '))} ? . anc ) w hen a measurement is 
performed, the probability of obtaining a measurement out¬ 
come must be independent of the eigenvectors not to break 
coherence between the eigenvectors. 

We divide the quantum circuit into seven steps from step (i) 
to step (vii) as shown in Fig. [13] We show the detail of how 
the eigenvectors are transformed after each step. 

First, we prepare a three-qubit input state 

|^ 0 ) )i, 3 |0>2 (104) 


in step (i), where we denote the index of the qubit corre¬ 
sponding to the first horizontal wire as 1 and that of the others 
likewise. After applying Hadamard gates in step (ii), we obtain 

ff 1 TT 3 |'J' (i) }i,3|+>2, (105) 

where |±) = -^=(|0) ± |1)). After applying Ci j3;2 in step (iii), 
we obtain 


1 

71 


(ff 1 /f 3 |^) 1 , 3 |o) 2 + z 1 fr 1 z 3 j? 3 |^) li3 | 1 ) 2 ). 


(106) 


After applying Hadamard gates and Pauli X operations in step 
(iv), we obtain 

—(XiX 3 |T' (i) )i, 3 |+)2 + I^M—>: 


(i) (ii) (iii) (iv) (v) (vi) (vii) 



Fig. 13. A protocol to implement a two-qubit unitary operation 
U g lobal( x 5 y 5 z ) over the butterfly network. We consider 7 steps presented 
in the quantum circuit and denote the steps by Roman numerals, (i) to (vii). 
The symbols of gates of the circuit are same as the ones given for Fig. [9] 


|^ W )i,3|0) 2 0 = 0,2) (107a) 

-|^ W )r,3|l)2 0 = 1,3). (107b) 


After applying C[ 3 ., in step (v), we obtain 

( e i(-y+*)|$(°)) li3 |o) 2 0 = 0) 

I ie < («'+*)|*( 1 )> li3 |l> 2 0 = D 

| e i <*'-*>|*< 2 >) 1 , 3 |0} 2 0 = 2) 

l *e < (-«'-*)|^( 3 )) li3 |l> 2 0 = 3). 


(108a) 

(108b) 

(108c) 

(108d) 


After applying a single qubit unitary operation u(x) given by 


uM = 2_0“ - ie ~“ ) 

' V 2 V / 

in step (vi), we obtain 

f A J l’f'«) 1 , 3 |+) 2 0 = 0,2) 

1 A J |T«) 1 , 3 |-) 2 0 = 1,3). 


009) 


(110a) 

(110b) 


After applying the projective measurement in the computa¬ 
tional basis and conditional unitary operations in step (vii), 
we obtain 

Ai |* w) ) ll3 dll) 


for any measurement outcome. We can verify that the prob¬ 
ability of obtaining a measurement outcome is A irrespective 
of eigenvectors. 


Appendix E 

Analysis of a bipartite property of four qubit 

STATES 

We prove that there is no pure state of four qubits |$)i, 2 ,3,4 
satisfying 

SCH#?$(|$» = 4, (112) 

SCH#^(|$» = 2, (113) 

SCH#^(|$)) = 2. (114) 

In I30l, it is shown that any pure states of four qubits can, up 

to permutations of the qubits, be transformed into one of the 
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following nine families of states by determinant 1 SLOCC: 

I*1> = ^(10000) + 11111)) + ^(|0011) + 11100)) 
+ ^( 10101 ) + 11010 )) + ^( 10110 ) + | 1001 )) 
l* 2 ) = ^(10000)+ 11111))+ ^(10011)+ 11100)) 

+c(|0101)+ 11010))+ |0110) 

|$3> = a(|0000) + 11111)) + 10101) + 11010)) 

+ | 0110 ) + | 0011 ) 

|$ 4 ) = a(|0000)+ 11111))+ ^-^(10101)+ 11010)) 

+~^~(|ono) + 11001)) 

+ -5=(|0001) + |0010) + |0111) + |1011)) 

v 2 

|$ 5 > = a(|0000) + |0101) + 11010) + 11111)) 

+i|0001) + |0110) — «|1011) 

|$ 6 ) = a(|0000) + |1111>) + |0011> + |0101> + |0110> 

|$ 7 ) = |0000) + |0101) + |1000) + |1110) 

|$ 8 ) = |0000) + |1011) + |1101) + |1110) 

|$ 9 ) = 10000) +10111), 

where a, b. c, d are complex parameters. 

Since the Schmidt number of a state cannot be in¬ 
creased under SLOCC and determinant 1 SLOCC is in¬ 
vertible, the Schmidt number of a state is invariant under 
determinant 1 SLOCC. Thus, we show that no state of the 
nine families simultaneously satisfies Eqs. (1112I) - (U 141) . There 
are three ways to divide four qubits into a pair of two 
qubits. We denote the set of Schmidt numbers of a four 
qubit state |4>) for all bipartite devisions as SCH#(|<f>)) = 
{SCH#^(|$)), SCH#^(|$)), SCH#^(|$))}. 

Theorem 5. There is no four qubit state +) £ 7Ti 0 TO 0 


7T.3 0 TO such that 

SCH#(|$)) = {4,2,2}. (115) 

Proof. By calculating the Schmidt rank for all bipartite devi¬ 
sions, we can easily check that 

Sch#(|$ 6 )) = {n 6 ,n 6 ,n 6 } (116) 

ScH#(|<f> 7 )) = {3,3,3} (117) 

Sch#(|$ 8 )) = {3,3,3} (118) 

Sch#(|$ 9 )) = {2,2,2}, (119) 


where n 6 = # jv^, + 4|a| 2 + |, \^1 + 4|a| 2 - 

and fS is the number of non-zero elements of set S. Since 
rig = 2 or ri6 = 3, these four states do not satisfy Eq. ( 11151 ). 

An element of SCH#(|<b 5 )) is # {l, y/2, 2|a|}. To satisfy 
Eq. (HH, a = 0 is required. Then 

Sch#(|$ 5 )) = {2,3,3}, (120) 

which does not satisfy Eg. (1115k 

An element of Sch#(|$ 4 )) is #{|&|} + #{x\x 3 — (3|a| 2 + 
2)x 2 + (3|a| 4 + 2|a| 2 + l)x — |a| 6 = 0}. To satisfy Eq. ( Il l 5b . 


the element must be 2 or 4. If the element is 2, since #{x\x 3 — 
(3|a | 2 + 2)x 2 + (3|a | 4 + 2|a | 2 + l)x — |a | 6 = 0} is larger than 
1 and is 2 if and only if a = 0 , we have 

a = & = 0 . ( 121 ) 

Then SCH#(|i> 4 )) = {2,2,2}. Thus, the element must be 4. 
Since jf{x\x 3 — (3|a | 2 + 2)a ; 2 + (3|a | 4 + 2|a | 2 + l):r— |a | 6 = 0} 
is 3 if and only if a 0, we have 

a^0,bf0. (122) 

Another element of SCH#(|<b 4 )) is #{\a — 6 |} + #{x\64x 3 + 
(• • • )a ; 2 + (• • • )x— |a— &| 4 |3a + &| 2 = 0}, where we abbreviate 
coefficients of X 2 and x. Since this element must be 2, it is 
necessary that 

a — b = 0 or 3a + b = 0. (123) 

The other element of Sch#(|$ 4 )) is #{\a + b\} + #{x\64x 3 + 
(• • • )x 2 + (• • • )x— |a + 6| 4 |3a— b\ 2 = 0}, where we abbreviate 
coefficients of x 2 and x. Since this element must be 2, it is 
necessary that 

a + b = 0 or 3a — b = 0. (124) 

We can easily check that it is impossible to simultaneously 
satisfy Eqs. < 1 1 22b - (l 1 24b . 

Sch#(|$ 3 )) is { 713 , 713 , 713 }, where 

713 = #{^2, |a + b\,\a - 6|}, (125) 

7/3 = #{Vl + 4|a| 2 ±l,Vl+4|6| 2 ±l}. (126) 

To satisfy Eq. (11 15k n' 3 must be 2, that is a = b = 0. Then 
713 = 1, which does not satisfy Eq. (1115k 
Sch#(|$ 2 )) is { 712 , 712 , 712 }, where 

n 2 = #{|a|,| 6 |,v/l + 4|c| 2 ±l}, (127) 

n' 2 = #{|a + &± 2 c|,Vl + |a- 6 | 2 ±l}, (128) 

ti 2 = #{|a - b ± 2c|, sjl + |a + 6| 2 ± 1}. (129) 

In the following, we verify that { 712 , 712 , 712 } cannot be 


{4,2,2}, {2,4,2} or {2,2,4}. 

1 ) {n 2 , n' 2 , Ji 2 } 7 ^ {4, 2, 2}: 

If n 2 = 4, it is necessary that 

a 0, bf= 0, c^0. (130) 

If 7i2 = 2, it is necessary that 

a — b = a + b + 2c = 0, (131) 

a — b = a + b — 2c = 0, (132) 

or a + b — 2c = a + b + 2c = 0. (133) 

If ri 2 = 2, it is necessary that 

a + b = a— b + 2c = 0, (134) 

a + b = a-b-2c = 0, (135) 

or a — b — 2c = a — b + 2c = 0. (136) 


We can easily check that it is impossible to simultane¬ 
ously satisfy Eqs. ( 1 1 30b -( l 1 36b . 

2 ) { 712 , 71 , 2 , 71 , 2 , }^{ 2 , 4 , 2 }: 
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If n 2 = 2, it is necessary that 

a = 6 = 0, (137) 

a = c = 0, (138) 

or 6 = c = 0. (139) 

With the necessary condition for n 2 = 2, we obtain that 

o = 6 = c = 0. (140) 

Then, it is impossible to satisfy n 2 = 4. 

3) {n 2 ,n' 2 ,n' 2 '}^ {2,2,4}: 

If n 2 = 2, it is necessary that 

o = 6 = 0, (141) 

o = c = 0, (142) 

or 6 = c = 0. (143) 

With the necessary condition for n' 2 = 2, we obtain that 

a = 6 = c = 0. (144) 

Then, it is impossible to satisfy n 2 = 4. 

Finally, we analyze Sch#(|$i)). SCH#(|$i)) is 

(n-i, n \, n"}, where 

«i = #{|a|,|6|,|c|,|d|}, (145) 

n'i = #={|a + 6 — c — d|, |a — 6 + c — d\, 

| — a + 6 + c — d|, |a + 6 + c + d|}, (146) 

n" = #{| — a + 6 + c +d|, |a — 6 + c + d|, 

|o + 6 — c + d|, |a + 6 + c — d|}. (147) 

Note that n\, n\ and n" are invariant under permutation of 
a, 6, c and d. We verify that {rti, n\ , n"} cannot be {4, 2, 2}, 
{2,4,2} or {2,2,4} in the following. 

1) {«!,<,<} ^{4,2,2}: 

If rti = 4, it is necessary that 

a / 0, 6^0, c^0, d^0. (148) 

If n'| = 2, it is necessary that in general 

a + 6 — c— d = 0, o — 6 + c — d — 0 (149) 

a = d, 6 = c. (150) 

Then 

n" = #{|26|,|2o|,|2o|,|26|} = 4. (151) 

2) {m,n},n"} ± {2,4,2} and {m,n},<} ± {2,2,4}: 

If ni = 2, it is necessary that in general 

o = 0, 6 = 0, c^O, d^O. (152) 

Then 

n[ = n" = #{|c + d|, |c + d|, |c - d|, |c - d|}. (153) 

□ 


Appendix F 

A CLUSTER NETWORK WITH LOOPS 
A cluster network with loops is defined as follows. 

Definition 4. A network G = {V, £,X, O} is a generalized 
cluster network if and only if for some k > 1 and N > 1, 

V = {t+j; 1 < * < fc, 1 < j < N} 

X = {i+i; 1 < i < k} 

O = {v itN -, 1 < i < k} 

E = S sub U/C (154) 

where 

Ssub + S corn p , 

Scomp — {(^mjl 1 5^ <7l<fc,l<J< A^)}, 

£ = {(fi+t+j+i); 1 < i < k, 1 < j < N- 1)}. 

(155) 

For this network, if there exists a loop of vertical edges 
£ C Ssuf, such that for some j, L and {i m }^ =1 , 

E = {ei = (vi 1 j,Vi 2 j),e2 = { v i 2 ,ji v i3,j')i 

j _j ) |c m e n if nr ^ ti}(156) 

it allows to perform a cyclic permutation that transmits a 
qubit state from Vi lt j to Vi 2 j , from Vi 2 j to ti,; 3 j and so on 
(by consuming Bell pairs corresponding to the looped vertical 
edges for teleportation), in addition to performing controlled 
unitary operations presented in Section IV. Thus, quantum 
computation over a cluster network with loops of vertical 
edges may have more capability than that without a loop. Note 
that a condition for the implementable unitary operations over 
this type of cluster networks with loops are still restricted by 
Theorem 2 and 4. An extension of our results for more general 
networks is an open problem. 
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